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d , Abstract. A recently developed algebraic approach for constructing coherent states 

for solvable potentials is used to obtain the displacement operator coherent state of the 
Poschl- Teller potential. We establish the connection between this and the annihilation 
operator coherent state and compare their properties. We study the details of the 
^vq ■ revival structure arising from different time scales underlying the quadratic energy 

J> , spectrum of this system. 
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1. Introduction 



Since its introduction by Schrodinger [TJ, coherent states (CSs) have attracted 
considerable attention in the literature [21 [31 HI El EJ [7J . A variety of coherent states, 
e.g., minimum uncertainty coherent state (MUCS), annihilation operator coherent state 
(AOCS), displacement operator coherent state (DOCS) and recently Klauder type CS 
[I], possessing temporal stability, have been constructed and applied to diverse physical 
phenomena [3j. Coherent states of systems possessing nonlinear energy spectra are of 
particular interest as their temporal evolution can lead to revival and fractional revival, 
leading to the formation of Schrodinger cat and cat-like states P, EJ [10] . A celebrated 
example in quantum optics of the aforementioned phenomenon is the coherent state 
in a Kerr type nonlinear medium In quantum mechanical potential problems, 

Hamiltonians for potentials like, Poschl- Teller, Morse and Rosen-Morse (RM) lead 
to nonlinear spectra. Time evolution of the CSs for these potentials, a subject of 
considerable current interest [121 H31 HH LTSl HE! [171 UHl LTE EQl EH E2], can produce 
the above type of states. 

The simplest way to construct CSs is a symmetry based approach [5J. It is well- 



ex, 



1, one can construct all the 



known that making use of the Heisenberg algebra 
above type of CSs for the Harmonic oscillator, which are identical to each other. In 
many physical problems, groups like SU(2) and SU(1,1) manifest naturally, enabling a 
straightforward construction of CSs. For the identification of the symmetry structure 
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of quantum mechanical potential problems, recourse has been taken to a number of 
approaches, starting from the factorization property of the corresponding differential 
equations. For the so called shape invariant potentials, super symmetric (SUSY) 
quantum mechanics [23] based raising and lowering operators have found significant 
application. A Klauder type CS, using a matrix realization of the ladder operators, has 
also been constructed @]. The fact that, SUSY ladder operators act on the Hilbert space 
of different Hamiltonians, has led to difficulties [15] in proper operator identification 
of the symmetry generators. The ladder operators have been taken to be functions of 
quantum numbers, which makes the corresponding algebraic structure ambiguous. This, 
in turn, creates difficulty in establishing a precise connection between the complete set 
of states describing the CS and the symmetry of the potential under consideration. 
In a number of approaches an additional angular variable have been employed [24] to 
identify SU(1,1) type algebras for describing the infinite number of states of some of 
these potentials. Taking advantage of the shape invariance property, quantum group 
type algebras have also been used for describing the Hilbert spaces [25J. 

Recently, a general procedure for constructing CSs for potential problems have 
been developed by some of the present authors [26]. The approach makes use of novel 
exponential forms of the solutions of the differential equations associated with these 
potentials for identifying the symmetry generators [27J. No additional variables are 
introduced and unlike SUSY based approaches one stays in the Hilbert space of a given 
quantum problem while unravelling its symmetry structure. The present paper makes 
use of this approach to study the DOCS of the Poschl- Teller potential. The primary 
motivation for considering the Poschl- Teller potential is two-fold. First of all, it has 
a quadratic spectrum leading to a rich revival structure for its CS, which can lead 
to the formation of cat-like states. Secondly, many other potentials can be obtained 
from the Poschl- Teller potential by appropriate limiting procedure and point canonical 
transformations. Hence, the CSs obtained in this case may have relevance to other 
potentials. The temporal evolution, auto-correlation and quantum carpet structures 
[HI Uni EH] of the CSs are carefully analyzed for delineating their structure and various 
time scales present in this problem. These properties are then contrasted with the 
corresponding ones of the AOCS [26] . 

The paper is organized as follows. In the following section, we briefly outline 
the procedure to identify the symmetry generators for quantum mechanical potential 
problems, based on hypergeometric and confluent hypergeometric equations. These 
symmetry generators are then used for constructing the DOCS for general quantum 
mechanical potential problems. Dual nature of the DOCS with the AOCS is algebraically 
established. In Section 3, the DOCS for the Poschl- Teller potential is constructed and 
its properties studied. We identify and analyze the various time scales of the system in 
Section 4 and compare the quantum evolution of the CS with the classical motion. We 
conclude in Section 5, after pointing out various directions for further work. 
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2. Algebraic structure of quantum mechanical potential 
problems 



As is well-known, the Schrddinger equation for a number of solvable potentials can 
be connected with the hypergeometric (HG) and confluent hypergeometric (CHG) 
differential equations (DEs). For example, harmonic oscillator, Coulomb and Morse 
potentials belong to the CHG class, whereas Poschl- Teller and Rosen-Morse belong 
to the HG class. Below, we briefly outline the steps of a novel procedure for solving 
DEs which connects the solution of a DE with the space of monomials [29]. This 
is subsequently used for identifying the symmetry generators underlying quantum 
mechanical potential problems. 

A single variable linear differential equation can be easily cast in the form, 

[F(D) + P(x, d/dx)]y(x) — (f) 

where the first part F(D) is a function of the Euler operator D = xd/dx, possibly 
including a constant term and P(x, d/dx) contains all other operators present in the DE 
under study. The solution can be written in the form, 

-P(x, d/dx) 



n=0 



(2) 



F(D) 

with the constraint F(D) x x = [29]. Using Eq. (T5]) the polynomial solutions of the 
HG and CHG can be written in closed form exponential forms 



and 



2 F l( -n, b; c; x) = (-1)" T f + n ™ e vh«''& a » , (3) 

r(c + n)T(b) 

xi^-n; c; x) = (-1)" -p^-e p ^x n . (4) 

1 (c + n) 

The exponential forms of these solutions are ideal for identifying algebraic structures of 
the solution spaces. For that purpose, one first identifies raising and lowering operators 
in the space of monomials. The operators at the level of polynomials can be obtained 
through similarity transformations. The simplest lowering operators at the level of 
monomial for CHG and HG functions can be taken (271 as 



K — x^ 2 +c d and K — 1 (x ^ + c ^ ) (5) 
dx 2 dx' (D + b) dx 2 dx ' 

respectively. The only criterion in choosing these operators at the monomial level is 
that, these do not lead to divergent expressions after the similarity transformation. It 
can be easily shown that, for the CHG case, the following generator form a SU(1,1) 
algebra at the monomial level: 

d 2 d d c , . 

A_ = x-~^r + c— , K + = x, andK 3 = x— + -. (b) 
dx 1 dx dx I 
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Similarly, for the HG case, the SU(1,1) generators are given as, 

- 1 , d 2 d 

~ (D + by x d^ + c di )} 

d c 

K , = (D + b- l)x, and K = x— + -. (7) 

dx 2 

Modulo a normalization, the DOCS for the HG type DE can be written, at the 
monomial level, as 

&( x ) = e P R +x°, (8) 

Here, x° = 1 is the fiducial state satisfying, 

K„x° = 0. (9) 

To find the CS x{ x i ft) a t the level of the polynomial, we make use of the exponential 
form of the solution in Eq. (j3j). The DOCS, x{ x ift)i can then be written as, 

X (x,ft) = e- R -e m +x° 

oo on 

= e~ K - ^r[(D + b-l)x] n x 



n=0 U - 



t P n . T{b J n) c- K -x" 



n=0 



n\ T{b) 



E^(-D -j^^-n^cx). 



(10) 



It is worth noting that, since the similarity transformation does not affect the algebraic 
structure, the SU(1,1) algebras remain intact at the polynomial level, albeit with 
different expressions for the generators. 

It is interesting to note that, at the monomial level, the DOCS found above is 
nothing but the AOCS of i.e. K-^{x) = ft^ ,3 (x), where 

~ 1 ( d?_ d_\ 

~ (D + b)(D + c) [ X dx 2+C dxr { } 



One notices that 



K_,K 4 



1. Hence, the above procedure is akin to the oscillator 



construction of AOCS. We can also identify a K + , which leads to the oscillator algebra 



1: 



(D + b-l\ , . 

k +={dT7^t) x - (12) 

The AOCS considered earlier [26J, is the eigen state of K_ and is of the form e /3K+ x°. 
This relationship between DOCS and AOCS has been referred earlier as duality of these 
two type of CSs [30] . Thus far, the specific nature of the potential has not been invoked. 
Now, we shall use this form to find out the CS for Poschl- Teller (PT) potentials. 
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3. Coherent state for the symmetric-Poschl-Teller potential 



The trigonometric Poschl- Teller potential belongs to the HG class having an infinite 
number of bound states. Hence it is natural to expect an underlying SU(1,1) algebra as 
its spectrum generating algebra. In reference [26] AOCS of the Poschl- Teller potential 
has been constructed, making use of a novel exponential form of the solution of the 
hypergeometric differential equation. Below we will concentrate on the construction of 
DOCS, following the same procedure and study its properties. We also compare the 
properties of DOCS and AOCS. 

The eigen values and eigen functions [31] of the symmetric-Poschl-Teller potential 



VspAv) 



2^,2 



h a 
2m 



P(P-1) 
cos 2 ay 



p > 1, 



(13) 



are given by, 



E 



SPT 



2^,2 



h a 



{n + pf 



n 



0,1,2, 



SPT, 



X) 



2m 

' a(n\)(n + p)T(p)T(2p) 
^T(p+\)T(n + 2p) 



and 



(l-x 2 )i^(x), 



with x = sin ay. Using the relation [32] 

V{2p + n) 



C>(1 - 2s) 



2-Fi(-n, b, c,x) 



and x 



T(2p)r(n + 1 

1 — 2x, where b = 2p + n and c = p + 1/2, we obtain from Eq. ([TO 



E 

n=0 



nl 



T(p + n + l/2)r(2p) 



C p n {x) 



(14) 



(15) 



(16) 



T(2p + n)T(p+ 1/2) _ 

Now multiplying Eq. ( TTBT) by (1 — x 2 ) p /2 and comparing with Eq. (fl4l) . we get the 
coherent state in energy eigenfunction basis as, 



where 



n=0 



SPT7- 



x) 



(17) 



r(p + n + l/2) S 



1/2 



_T(2p + n)T(n + l)(n + p)_ 
For comparison, the eigen function distribution for AOCS can be written as 

1/2 



d 



AOCS 



(7)" 



1 



(19) 



_r(2p + n)T(n + l)(n + p)_ 

We can also obtain the DOCS of a general trigonometric Poschl- Teller potential [1] 
modulo a normalization factor, in the same manner as for the symmetric-Poschl-Teller 

case : 

oo 

/>T (^) = E d p J^ T (x) (20) 

n=0 



X 
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Figure 1. The \d n \ 2 plots of DOCS and AOCS of symmetric-Poschl- Teller potential 
for p = 15. 



where 



PT 



d 



T(k + n + l/2)T(p + n + l/2) 
(k + p + 2n)T(n + l)T(k + p + n) _ 

1 1/2 



1/2 



2a(k + p + 2n)T(n + l)T(k + p + n) 
T(k + n + l/2)T(p + n + l/2) 

x (1 - x) p/2 (x) k/2 Pt 1/2 ' p ~ 1/2 (l ~ 2x). 



(21) 



Although the symmetric-Poschl- Teller potential has infinite number of bound 
states, only a few states contribute appreciably to the sum, which is peaked around 
n = n. In Fig. [TJ we compare the nature of the distributions of the eigen states 
for AOCS and DOCS. For the purpose of comparison, we have taken the coherence 
parameters such that, the distributions are comparable. It is found that, both the eigen 
state distributions, peaked at n = 9, involve the same eigen states (from n=0 to n=30) 
for the same potential (p = 15). For AOCS, the distribution resembles a Gaussian 
distribution and is more sharply peaked, as compared to the DOCS. Larger (3 value 
makes the distribution flatter for DOCS, as seen in the dashed curve in Fig. [lj We now 
proceed to study the spatio-temporal dynamics of these wave packets. 

4. Revival dynamics of coherent state 

The time evolution of CS x{%ifl) can be written as 

oo 

X {x,t) = J2 d nMx)e~ iEnt - (22) 

n=0 
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Figure 2. Probabity density plot of DOCS of symmetric-Poschl- Teller potential at 
different times for (3 = 0.8 and p = 10. 



As the energy expression contains terms up to n 2 , the system shows revival and fractional 
revival but no super-revival phenomenon. All graphs are plotted in time, scaled by the 
revival time T rev = An /a 2 , . 

In order to throw more light on the structure of the revival pattern, we note that 
the eigen functions satisfy 

ip n {-x) = {-l) n ip n {x). (23) 

From Eq. (T2"2"j) . we can easily obtain the CS wave packet at time t = 1/2 T rev as 



|xO, t = -T rev )\ 2 = \x(-x,t = 0)| 2 . 



(24) 



Thus, at time t = 1/2 T rev , a mirror image of the initial wave packet is produced at the 
opposite end of the potential well (Fig. [2] and [3D- This can be observed as a bright ray 
at time t = 1/2 T rev , in the quantum carpet structure (Fig. H]). The auto-correlation 
function 



yields 



A(t) = {x{x,t)\x{x,0)), 



A(t = 1/2) = £ \d n \ 2 - K\ 2 

(A+n) even (A+n) odd 



(25) 
(26) 
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Figure 3. Probabity density plot for AOCS of symmetric-Poschl- Teller potential at 
different times for f3 = 30 and p = 10. 




0.5 



Figure 4. (Colour Online) Quantum carpet of the displacement operator coherent 
state of symmetric-Poschl- Teller potential for (3 — 0.8 and p — 10, brightness signifies 
the maximum. 



As d n oscillates rapidly, it will not contribute significantly to the |^4(t)| 2 (Fig. [5]), at 
t = 1/2 T rev . At time t = 1/4 T rev the CS wave packet becomes 



V2 
1 



-in/4 



xO,o) 



-x,0) 



\ X {x, 0)| 2 + |x(-x,0)p 



(27) 



In this case, the wave packet breaks up into two parts which are situated at the 
two opposite corners of the potential well (Fig. [2] and [3]). This gives rise to two bright 
spots at the two vertical ends of the quantum carpet at t — .25. At the same instance, 
the auto-correlation function only gives a peak, as manifested in Fig. [5j To explain the 
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Figure 5. Plot of square modulus \spt(x( x i t)\x( x > ®-))spt\ 2 of the auto correlation 
function of DOCS as a function of time, for (3 — 0.8, p = 10. 



probability density plot at t = 0.125, we consider a fictitious classical wave packet 



Xd(x,t) = ^2 d n ip n (x)e 



n=0 



(28) 



2n 



which behaves like the initial wave packet at small time (order of T c i) where T c i — 
Using the discrete Fourier transform (DFT), the original CS wave packet (Eq. ( 1221) ) at 
time t = ^Trev can be written as a linear combination of classical wave packet of Eq. (1281) 

as 

i-i 

r n 

(29) 



S p=0 s I 



where 



i-i 



a p = - ex P 

1 n=0 



2rri(n- — n 2 -) 



(30) 



I s- 

Here, r and s are two mutually prime integers and I is the period of the quadratic 
phase term. In general, I — |, if s is integral multiple of 4 and I — s in all other cases. 

s 
2 



In this case - = I and I — I- = 4, Substituting in Eq. (1291 and Eq. fl30|) . we get 



(0) (1) (2) 

«o Xci + a iXd' + a 2Xd 



(3) 
«3Ad 



where 



Xcl 



1 i 

Xc/(x, -T re „ + -T d ), i = 0,l,2,3 



(31) 
(32) 



and a = -a 2 = ^(1 - i) ; a x = a 3 = 1/2. 

The above expression (I3"T|) signifies that the wave packet has broken into four parts, 
each of them differing by a phase 7r/4. In the probability density 



L \ 2i I i. I _l I .1 I i I 

Xdl + \Xd\ + \Xd\ + \x. 



,2 1 2 



3 |2 



,4 |2 
dl 



+ iRe 



A^Xd^ 4 - Xd *xi + x i :xy^ 



2* 3 -i ff /4 , 2* 4 _ 3* 4 «r/4 
Ad Ad e + Ad Ad Ad Ad e 



(33) 



(34) 
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Figure 6. Plots of expectation values of x for (a) the coherent state of the Poschl- 
Teller potential, with different values of f3 where p = 5, k = 5 and a = 2.(b) classical 
solution of Poschl- Teller potential, with a — 0.25, m = 1, Vb = 4, k = 5, p = 5 and 
E c = average value of energy for (3 = 0.1. 



the first term carries the contribution from the individual subsidiary waves and the 
second term arises due to the interference between them. We note that the xli an d Xd 
are spatially well separated, giving very less contribution in interference. The dominant 
interference term is Xd Xd-> as Xd an d Xd are n °t spatially separated. Thus, at time 
t = 0.125, wave function splits into four parts, but their interference at the middle 
gives a strong peak, rather than giving four distinct similar waves. For comparison, the 
wave packet structure for AOCS of symmetric-Poschl- Teller potential is also shown in 
Fig. In this the initial wave packet is sharper, the interference terms are less 

dominant than that of DOCS. 

We have observed that the initial wave packet remains in the left corner of the 
potential well and oscillates due to the impulse from the well and at later times, as 
it spreads, being away from the boundary of the well. This is quite transparent from 
the quantum carpet structure, which gives the space time rays of probability density of 
the corresponding coherent states. We note that, the rays in quantum carpet are not 
straight lines which is the case for infinite square-well. 

In order to contrast the temporal evolution of the DOCS with the classical motion, 
we note that for a particle of energy E, 



x{t) = a arccos 



hvAcosi 

2 x m a 




(35) 



where A = (l-I( v ^7+ v ^)2)(i_I( v ^7_ v ^)2) 5 Vq = £ and ai = gp(p-l), fa = 
j^-k(k — 1); with the condition 

£ c >^(vW-l) + V^-l)) 2 - (36) 

This classical trajectory is shown in Fig. E](b). The expectation value of the position 
with respect to the coherent state of general trigonometric Poschl- Teller potential is 
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(37) 



where 

having 
A n = - 
and 

C n = (P) 2n 



N J2 ( 2A n cos [2 a 2 (2n + p + k + l)t] - C n 



n=0 



2n+l 



2T(p + n + 3/2)r(fc + n + 3/2) 



T(p+k + n)T(n + l)(2n + p + k)(2n + p + k + l)(2n + p + k + 2) 
T(p + n + l/2)r(fc + n + 1/2) (A: + p - l)(fc - p) 



r(p + fc + n)r(n + l)(2n + p + k)(2n + p + fc - l)(2n + p + k + 1) 

N being the normalization constant. This (x(t)) is plotted in Fig.[6](a) which nearly 
matches with the classical trajectory for very small values of (3 (solid line in Fig. E](a)). 
In this case only a few eigen states contribute to the coherent state wave packet. Sudden 
changes in the (x(t)) values are the signatures of revivals and the fractional revivals [33 J. 



5. Conclusions 



In conclusion, the algebraic procedure used here, for constructing CSs for potentials 
based on confluent hypergeometric and hypergeometric differential equations depends 
on the fact that, the solutions of the above differential equations can be precisely 
connected with the space of monomials. This leads to a straightforward identification 
of symmetry generators, without taking recourse to additional angular variables or 
SUSY type multiple, related Hamiltonians. The nature of the specific potential enters 
through the corresponding ground states and by fixing the parameters and variables of 
the above solutions. We have concentrated here on the CS of Poschl- Teller potential, 
since various potentials can be obtained from the same, through limiting of parameters 
and point canonical transformation. The time evolution of the CS for this potential, 
having non-linear spectra, produces cat like states in fractional revivals. We contrasted 
the properties of the two different types of CSs possible here, as well as the temporal 
evolution of the CS, with classical motion. As has been noted earlier, this procedure 
easily extends to more complicated non-linear coherent states arising from deformed 
algebras, a subject we intend to take up in future. We also would like to analyze the 
subject of mesoscopic superposition and sub-Planck scale structure [33], possible in this 
type of quantum systems. 
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